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■ Cellular membranes are a heterogeneous mix of lipids, proteins 
, and small molecules. Special groupings of saturated lipids and 

■ cholesterol form a liquid-ordered phase, known as 'lipid rafts,' 
serving as platforms for signaling, trafficking and material transport 
throughout the secretory pathway. Questions remain as to how 
the cell maintains heterogeneity of a fluid membrane with multiple 

■ phases, through time, on a length-scale consistent with the fact 
, that no large-scale phase separation is observed. We have utilized 

■ a combination of mechanical modeling and in vitro experiments to 
show that membrane morphology can be a key player in maintaining 

" this heterogeneity and organizing such domains in the membrane. 
We demonstrate that lipid domains can adopt a flat or dimpled 
morphology, where the latter facilitates a repulsive interaction that 
slows coalescence and tends to organize domains. These forces, that 

■ depend on domain morphology, play an important role in regulating 
lipid domain size and in the lateral organization of lipids in the 

■ membrane. 

The plasma and organelle membranes of cells are composed 
of a host of different lipids, lipophilic molecules and membrane 
proteins Together, they form a heterogeneous layer ca- 
pable of regulating the flow of materials and signals into and 
] out of the cell. Lipid structure and sterol content play a key 
role in membrane organization, where steric interactions and 
energetically costly mismatch in the hydrophobic structure of 
lipid tails result in lateral phase-separation. Saturated lipids 
and cholesterol are sequestered into liquid-ordered (Lo) do- 
mains, often known as 'lipid rafts', from an unsaturated liquid- 
disordered (L^) phase [U |3l 2]. Domains composed of sat- 
urated sphingolipids and cholesterol, with sizes in the range 
of 50 — 500 nm, have been implicated in a range of bio- 
logical processes from lateral protein organization and virus 
uptake to signaling and plasma-membrane tension regulation 
[SliElIHliaiinilllllllElIIlIISlIIlIIT] how the cell main- 
tains the lateral heterogeneity of lipids over time, and what 
physical mechanism might be responsible for the spatial orga- 
nization of these domains, challenges classical theories of phase- 
separation and 'domain ripening' (such as Cahn-Hilliard kinet- 
ics m]). The maintenance of lateral heterogeneity is thought 
to arise from a combination of lipid recycling and energetic bar- 
riers to domain coalescence |191I201[2T| (potentially provided by 
transmembrane proteins [22^), resulting in a stable distribution 
of domain sizes. The precise origin of this energy barrier and 
the nature of its dependence on membrane elastic properties 
remains unclear. 

The simplest physical model that describes the evolution of 
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lipid domain size and position predicts that domains diflFuse 
and coalesce, such that the number of domains constantly de- 
creases, while the average domain size constantly increases |18) . 
Indeed, models of two-dimensional phase separation have been 
studied in detail for many physical systems [231 [2^ 1251 [26] , 
and where the phase boundary is unfavorable and character- 
ized by an energy per unit length [27], the domain size grows 
continuously (oc t^^^) [181 [28l [29] . However, membranes can 
adopt three-dimensional morphologies that affect the kinetics 
of phase separation [30l [3L 32 , 33, jSl]. In those cases where 
morphology is considered as part of the phase separation model, 
novel coalescence kinetics emerge [31]. Experimentally, model 
membranes have shown that nearly complete phase separation 
on the surface of a cell-sized vesicle can be reached in as little 
as one minute This seems inconsistent with the fact that 
on the cell surface, much smaller domains persist on that same 
time-scale [21] and no large-scale phase separation is observed. 
With these facts in mind, our central questions are: how can 
membranes that have phase-separated maintain their lateral 
heterogeneity on long time scales and short length scales? Are 
there membrane-mediated (i.e. elastic) forces that inhibit co- 
alescence and spatially organize domains? 

We begin to answer these questions by examining the ener- 
getics of the membrane using a linear elastic model. A phase- 
separated membrane is endowed with bending stiffness, mem- 
brane tension, an energetic cost at the phase boundary, and 
domains of a particular size. Membrane bending and tension 
establish a natural length-scale over which a morphological in- 
stability develops that switches domains from a flat to 'dim- 
pled' shape, similar to classical Euler buckling [35] (see Figure 
[1]). The dimpling instability is size-selective and 'turns on' 
a membrane-mediated interaction that inhibits domain coales- 
cence. This transition is a precursor to budding, and is distinct 
from transitions that require spontaneous curvature. While 
variations in membrane composition may change specific pa- 
rameter values, the mechanical effects we describe are generic. 
Thus, these systems exhibit shape-dependent coarsening kinet- 
ics, that are relevant for a broad class of two-dimensional binary 
fluid systems. The interaction between domains is a mechanical 
effect, and we use a model treating dimpled domains as curved 
rigid inclusions to distill the main principles governing this in- 
teraction. The confluence of membrane properties required for 
this morphological change and its attendant forces lies squarely 
in the biological regime. Experimentally, we use a model mix- 
ture of lipids and cholesterol to show that such an interaction 
exists between dimpled domains and is well approximated by 
a simple model. Combined with lipid recycling [19], we of- 
fer elastic interactions as a mechanism for the maintenance of 
lipid lateral-heterogeneity and organization of domains in cel- 
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Fi glire ll Three dimensional rendering of a dimpled lipid domain 
in dimensionless coordinates. For a domain (shown in red), a com- 
petition between bending, membrane tension, and phase boundary 
line tension results in a morphological transition from a flat to a dim- 
pled state as depicted above. The dimple costs bending energy but 
relieves line tension by reducing the phase boundary length (shown 
as a white line around the domain). This morphology facilitates in- 
teractions between domains that signiflcantly alter the kinetics of 
coalescence and lateral lipid organization. The projected domain 
radius is po = To 1)^2 • 



lular membranes. 

The first section of the paper outlines the energetic contri- 
butions to the mechanical model, and predicts the conditions 
under which domain dimpling occurs. The second section 
outlines how dimpled domains facilitate an elastic interaction 
and compares the model interaction to our measurements 
made in phase-separated giant unilamellar vesicles. 

The Elastic Model and Morphological Transitions 

The energetics of a lipid domain are dominated by a competi- 
tion — on one hand the applied membrane tension and bending 
stiffness both energetically favor a flat domain; on the other 
hand the phase boundary line tension prefers any domain mor- 
phology (in 3D) that reduces the boundary length. We use a 
continuum mechanical model that couples these effects, relating 
the energetics of membrane deformation to domain morphol- 
ogy. As we will show, this competition results in a morpholog- 
ical transition from a flat to dimpled domain shape, where two 
dimpled domains are then capable of interacting elastically. 

Lipid domains in a liquid state naturally adopt a circular 
shape to minimize the phase boundary length [2 , allowing us 
to formulate our continuum mechanical model in polar coor- 
dinates. We employ a Monge representation, where the mem- 
brane mid-plane is described by a height function h{r) in the 
limit of small membrane deformations (i.e. |V/i| < 1). With 
this height function, we characterize how membrane tension, 
bending, spontaneous curvature and line tension all contribute 
to domain energetics. 

Changes in membrane height alter the projected area of the 
membrane and hence do work against the applied membrane 
tension, resulting in an increase in energy written as 

Ctens = TTT (^J (V/ii)Vdr + J (V/i2)VdrJ , (1) 

where r is the constant membrane tension, ro is the projected 
radius of the domain, hi is the height function of the domain 
and /i2 is the height function of the surrounding membrane [361 
I37j . Membrane curvature is penalized by the bending stiffness 
with a bending energy written as [39l [36] 

Gbend =7r4^^ (^a J^^ {y^hi)\dr + (V'/is)' rdr^ . (2) 



Our model allows the domain and surrounding membrane to 
have differing stiffnesses, and ^t^^^ respectively, character- 
ized by the parameter a = ^5^V^^^ from this point on 
we drop the superscript on k,^^\ Recent experiments suggest 
that the bending moduli of a cholesterol-rich domain and the 
surrounding membrane are roughly equal [H [38] , and hence for 
simplicity, we assume the bending moduli of the two regions 
are equal {i.e. a = 1), unless otherwise noted. In addition to 
bending stiffness, the domain may exhibit a preferred 'spon- 
taneous' curvature due to lipid asymmetry or protein binding 
[401 134] . The contribution of domain spontaneous curvature 
can be simplified to a boundary integral, which couples to the 
overall curvature field by 

pro 

Cspont = -27iaKbCo / (V^/ii) rdr = -27iaKbCoroe, (3) 
Jo 

where Co is the spontaneous curvature of the domain and e is 
the membrane slope at the phase boundary as shown in Figure 
[T] Further, we assume the saddle-splay curvature moduli are 
equal in the two regions, yielding no dependence on Gaussian 
curvature. In principle, this contribution could be accounted 
for with a boundary term, explored in detail in the supple- 
mentary information (SI). The phase boundary line tension is 
applied to the projected circumference of the domain, as shown 
in Figure [1] by Gune = 27rro7 where 7 is the energy per unit 
length at the phase boundary. 

Finally, a constraint must be imposed that relates the ac- 
tual domain area. A, to the projected domain radius ro- The 
energetic cost to change the area per lipid molecule is high 
(- 50 - lOOkBT/nm^ where ks = 1.38 x 10"^^ J/K and 
T = 300 K [U), hence we assume the domain area is conserved 
during any morphological change (see SI for details). We im- 
pose this constraint using a Lagrange multiplier. To, with units 
of tension by 

Carea = To (^TT ^ (V/ll) Vdr + TTr^ - . (4) 

This results in an effective membrane tension within the do- 
main Ti = T + To, which must be negative to induce dimpling. 
Examining the interplay between bending and membrane ten- 
sion, we see that two natural length scales emerge - within the 
domain we define Ai = GKhjrx and outside the domain we 
define A2 = \/ t^h/r. These length scales allow us to define the 
relevant dimensionless parameters in this system. 

The total free energy of an elastic domain and its sur- 
rounding membrane is then the sum of these five terms, G — 

Gtens + Ghend + ^spont + Gune + Carea ■ Dctalls OU all thc tCrmS 

in the free energy can be found in the SI. With this free energy 
in hand, we examine how the morphology of a circular domain 
evolves as we tune domain size and the elastic properties of the 
membrane. 

The height field and radius can be rescaled by the elastic 
decay lengths such that the Euler- Lagrange equation for the 
domain can be written in the parameter- free form V^(V^ + 
/3^)?7i = 0, while the equation for the surrounding membrane is 
V^(V^ — 1)?72 = 0, where the dimensionless variables are defined 
by X2r]i — hi, X2P = r, X2P0 — To and /3 = 1X2/ \i. Using the 
same dimensionless notation, the energy from line tension and 
spontaneous curvature can be written as Cune = 27ThibpoX ^-nd 
G^spont = -27TaK,bepoVo, with Vo = X2C0 and x = lX2/nh- The 
dimensionless line tension, x, is simply a rescaled version of the 
line tension 7 and is one of two key parameters that characterize 
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the morphological transition; the dimensionless domain area, 
a = ^/A|, is the second key parameter. 

The admissible solutions for 771 (p) and 772 (p) are zeroth or- 
der Bessel functions Jo{Pp) and i^o(p), respectively, with the 
boundary conditions |V?7i(0)| = |V?72(oo)| = and |V?7i(po)| = 
|V?72(po)| = e- The boundary slope, e, is the parameter that 
indicates the morphology of the domain; e = indicates a flat 
domain, while < |e| < 1 indicates a dimpled domain. The 
five contributions to membrane deformation energy yield a rel- 
atively simple expression for the total free energy, given by 
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Mechanical equilibrium is enforced by rendering the energy sta- 
tionary with respect to unknown parameters e, po, and (3^ 



^ = 0, ^ = 0, ^ = 0. 

de ' dpo ' df3 



(6) 



These equilibrium equations physically correspond to torque 
balance at the phase boundary, lateral force balance at the 
phase boundary and domain area conservation, respectively. 

Analysis of the equilibrium equations reveals a second- order 
transition at a critical line-tension, Xc, as shown in Figure (2] 
For X l^ss than this critical value, only the flat, trivial solu- 
tion with e = exists. At Xc a non-trivial solution describing 
buckled or dimpled morphologies emerges. For zero sponta- 
neous curvature, the bifurcation is defined by a transcendental 
characteristic equation 
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with P = (Xc/ po — ^)/cr and po — \foi~pK . For a given dimen- 
sionless domain area, a, this defines the critical line tension re- 
quired to dimple the domain. In Figure [2^ (inset), this relation 
is used to generate a morphological phase diagram that shows 
where in the space of dimensionless domain area and line ten- 
sion we find the discontinuous transition (i.e. bifurcation) from 
a flat domain, to a dimpled domain. Near the morphological 
transition the boundary slope grows as |e| oc \JxlXc — 1, indi- 
cating that a dimple rapidly deviates from the flat state. The 
transition is symmetric, in that both possible dimple curvatures 
have the same energy, and hence the domain is equally likely 
to dimple upwards or downwards. In the experimentally rele- 
vant limit of small dimensionless domain area, the complexity 
of eqn. [Tlis reduced to 



(8) 



This leads to the conclusion that the dominant parameter gov- 
erning domain dimpling at zero spontaneous curvature is X\^- 
For a small domain, the dimpling transition is directly reg- 
ulated by domain area, the bending modulus, and line ten- 
sion, but only weakly depends on applied membrane tension. 
Intuitively, domains dimple when line tension or domain size 
increase (subject to small a), as shown in Figure [2^(inset). 
Likewise, a decrease in bending stiffness, due, for instance, to 
changes in membrane sterol content [421 H3] , can also induce 
dimpling. The effects of applied membrane tension are weak 
because the change in projected area upon dimpling does not 
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Figure 2: Bifurcation diagram for dimpling transition at constant 
area (a = 7r/4, hit, = ^bksT, A2 = 500 nm, a = 1). Constant line 
tension and increasing area produces a qualitatively similar graph, 
a) At zero spontaneous curvature (vq = ^ black) the bifurcation 
is symmetric, the upper and lower branches are at the same en- 
ergy, and e = becomes unstable above the critical point (horizon- 
tal black dashed line). With finite spontaneous curvature (vo = 2, 
Co = (250 nm)"-*^ blue) the lower energy branch (upper) has non- 
zero e for all line tensions, asymptoting to the Vq = branch. At a 
line tension slightly higher than Xc for the Vo = case, a bifurca- 
tion yields a higher energy dimple with the opposite curvature as Vo 
(indicated by the second vertical dashed line). Inset: Equilibrium 
phase diagrams for a = 0.5 (red), a = 1 (green), and a = 2 (blue) 
(the dashed lines are the approximation of eqn. [8]) showing flat (F) 
and dimpled (D) domains, b) Energy difference between the flat 
and dimpled state, normalized by for domains with and without 
spontaneous curvature (vo = ^ black; Vq = 2 ^ blue). 



lead to a significant energy cost relative to the cost of bending 
and line tension. 

If membrane elastic properties are fixed (i.e. fixed Kb, r 
and 7), the dimpling-induced interactions 'turn on' only after 
a critical domain size is achieved. This scenario is encountered 
when two domains, too small to dimple on their own, diffusively 
coalesce into a larger domain capable of dimpling and hence in- 
teracting. Indeed, such a size-selective coalescence mechanism 
was observed recently in model membrane vesicles This 
constitutes a distinct class of coarsening dynamics, where clas- 
sical diffusion-limited kinetics are obeyed until the domain size 
distribution has matured past the critical size for dimpling - 
then domain coalescence is a relatively slow, interaction- limited 
process. 

For the model domain considered in Figure [2l with area 
a = 7r/4 (ro — 250 nm), the critical dimensionless line ten- 
sion is Xc — 13, corresponding to a critical line tension of 
7c 0.65 ksT /nm. This value compares well with both theo- 
retical estimates of the line tension |27[|45) . and the higher side 
of experimentally measured values [T, 46' Tf^ . 

Spontaneous curvature does not affect the Euler- Lagrange 
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Figure 3: Theoretical and experimental dimpled domain shapes. 
Domains are shown in red, surrounding membrane in blue, a) Min- 
imum energy dimples with and without spontaneous curvature (see 
legend, o; = 7r/4, cr = 1). b) Epi-fluorescence cross-section of a dim- 
ple on the surface of a GUV; the red and blue lines are a guide to the 
eye. c) ID model of interaction - dimples maintain shape, but tilt 
{(j)) as a function of separation distance {d). Dimples with the same 
sign of curvature repel, while dimples with opposite sign attract, d) 
Epi-fluorescence cross-section of two dimpled domains interacting on 
the surface of a GUV. Scale bars are 3 /xm. 



equations, and hence will not effect the class of equilibrium 
membrane shapes. However, domains with zero and nonzero 
spontaneous curvature exhibit qualitatively different behavior. 
Biological membranes can be asymmetric with respect to leaflet 
composition [421 [5] |48], endowing a domain with potentially 
large spontaneous curvature. The energetic contribution from 
spontaneous curvature takes the form of an additional line 
tension depending linearly on the slope taken by the domain 
boundary, e. This breaks the symmetry of the membrane, giv- 
ing an energetic preference to a dimple with the same curvature 
as the spontaneous curvature, and eliminating the trivial e = 
solution even at small line-tensions. As line tension increases, 
a bifurcation produces a second, stable, higher-energy dimple 
of the opposite curvature as Vq. The more energetically sta- 
ble branch of this transition corresponds to a dimpled state for 
all values of line tension and non-zero values of domain area, 
as demonstrated in Figure [2^. This predicts that as soon as a 
domain with finite spontaneous curvature forms, it dimples, re- 
gardless of size, and begins to experience interactions with any 
nearby dimpled domains. It is reasonable to expect that do- 
mains with similar composition will have similar spontaneous 
curvature, and hence form dimples whose curvature has the 
same sign. As we will show, dimples whose curvature has the 
same sign tend to interact repulsively. Such a mechanism of 
coalescence inhibition was observed recently in simulation [3l] . 

This indicates that control of spontaneous curvature via do- 
main composition or protein binding can regulate dimpling and 
hence domain interaction [491 148) . Indeed, recent experimental 
[50] and theoretical [51] work shows that protein binding and 
lipid asymmetry, respectively, lead to precisely these kinds of 
dimpled domains. 

Calculated shapes of dimpled domains induced by line 
tension and spontaneous curvature are shown in Figure [3^, 
alongside dimpled domains observed on giant unilamellar 
vesicles, shown in Figure [Sj) and d. 

Elastic Interactions of Dimpled Domains 

Given two domains that have met the criteria for dimpling, the 
deformation in the membrane surrounding the domains medi- 



ates an elastic interaction when they are within a few elastic 
decay lengths (A2) of each other. This equips us to begin ad- 
dressing how short length-scale and long time-scale membrane 
heterogeneity might be achieved. As previously stated, free 
diffusion sets the maximum rate at which a quenched mem- 
brane can evolve into a fully phase-separated membrane [18], 
where this evolution can happen in as little as a minute on 
the surface of a cell-sized vesicle [2]. On the other hand, re- 
cycling and hence homogenization of cellular membrane is a 
process that takes place on the time-scale of an hour or more 
[52) . Our measurements of domain interactions (detailed below 
and other data shown in SI) estimate the coalescence barrier 
between dimpled domains at 5 — lOksT. Hence, given the 
diffusion-limited rate of phase separation, interactions slow this 
process by approximately e"^ :^ 0.007 to e"^° ^ 0.00005. This 
makes the time-scale of lipid heterogeneity comparable to the 
time-scale of membrane recycling and even eukaryotic cell di- 
vision. 

The physical origin of domain interaction is explained by a 
simple model based on the assumption that the dimpled do- 
main shape is constant during interaction, but the domains are 
free to tilt by an angle 0, as shown in Figure This as- 
sumption was, in part, inspired by experimental observations 
of domain shapes on the surface of giant unilamellar vesicles, 
as shown, for example, in Figure [Sjl. The interaction energy 
is roughly an order of magnitude less than the free energy as- 
sociated with the morphological transition itself (see Figure 
[2)3), thus interaction does not perturb the domain shape sig- 
nificantly. Only allowing domains to rotate simplifies the inter- 
action between two domains to a change in the boundary con- 
ditions in the three regions of interest, shown in blue in Figure 
[3]:. Applying the small gradient approximation, the boundary 
slope is given by |e — 0| in the outer regions and by |e + (/)| in the 
inner region. With the single domain boundary slope, e, set by 
the energy minimization of the previous section {i.e. eqn. [6]), 
the pairwise energy is minimized at every domain spacing, 
by dG/d(j) = to find the domain tilt angle that minimizes 
the deformation energy (see SI for details). This results in two 
qualitatively distinct scenarios: two domains whose curvatures 
have the same sign repel each other, while two domains whose 
curvatures have the opposite sign attract each other. ScaUng 
arguments can be used to show that the strength of interac- 
tion between two dimpled domains increases roughly linearly 
with their area, so long as they are both larger than some criti- 
cal area (see SI for details) . Mathematically, the assumption of 
rigidly rotating dimpled domains on a membrane is identical to 
a previous 2D model of bending-mediated interactions between 
intramembrane proteins represented by rigid conical inclusions 

Independent of the effects of spontaneous curvature, slight 
osmolar imbalances and constriction due to the lipid phase 
boundaries create small pressure gradients across the mem- 
brane that tend to orient all dimples in a cell or vesicle in 
the same direction, resulting in net repulsive interactions be- 
tween all domains. Transitions between 'upward' and 'down- 
ward' dimples are infrequent, due to a large energy barrier. In 
the simplest case, where the domains are the same size, the 
tilt angle (j) monotonically increases as two domains get closer, 
(t){d) — ee~^. Likewise, the interaction energy, Knt(<^) — 
pi^~^ ^ increases monotonically with decreasing separa- 
tion. For direct comparison, we fit both the ID model outlined 
here and the 2D inclusion model [SSJ to the data of Fig(4l show- 
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^^^^^^^1 Center-to-Center Distance (fim) 

Figure 4: Measuring domain interactions on the surface of a vesi- 
cle, a) Three images of dilute interacting domains on the surface 
of the same vesicle (scale bar is 10 /im). b) The repulsive interac- 
tion potential of domains on the surface of the same vesicle as (a). 
The energy is measured in /c^T and distance is domain center-to- 
center. The blue dashed line is a fit to the ID interaction model in 
this paper, Vint{r) = aie ^/ 2 + (12, with elastic decay length 
A'^-^' = 240 nm. The orange dashed line is a fit to the model, 
VintW = 27r/t5 [(aia2)2Xo(r/A^^'^^) + a2a|K|(r/A^^'^^)] +04, with 

elastic decay length X^2^^ — ^70 nm, based on the theory of Weikl et 
al [53) . Both elastic decay lengths indicate a membrane tension of 
~ 4 X 10~^ /c^T/nm^. Errors bars are shown in green on the x-axis. 



ing that they are experimentally indistinguishable, though with 
a slightly different elastic decay length. 

To quantitatively compare our interaction model with exper- 
iment, we examined the thermal motion of small domains on 
the surface of giant unilamellar vesicles, as described in 'Mate- 
rials and Methods.' Membrane tension was regulated by bal- 
ancing the internal and external osmolarity, giving us coarse 
control over the elastic decay length A2. Through time, the 
distance between every domain pair was measured and the 
net results were used to construct a histogram. The poten- 
tial of mean force as a function of distance between domains is 
shown in Figure HJ). We selected vesicles that had a low den- 
sity of approximately equal-sized domains, and thus generally 
the interactions were described by a repulsive pairwise poten- 
tial. Though areal density of domains and generic data quality 
varied in our experiments (see SI), all data sets exhibit the re- 
pulsive core of the elastic interaction. Multi-body interactions 
occur, though infrequently; their effect can be seen as a small 
variation in the baseline of Figure UJ), which is not captured 
by the pairwise interaction model. At high membrane tension, 
when we would not expect dimpled domains, we qualitatively 
verified that domains coalesce in a rapid manner as compared 
to our low tension experiments (data not shown). Other recent 
experiments have also observed repulsive interactions between 
domains on low membrane tension vesicles and the lack thereof 
on taut vesicles [44] . 

Our measurement of the pairwise potential allows us to es- 
timate elastic properties of the membrane. The elastic decay 
length was fit to the ID and 2D interactions models described 
above, and found to be X^2^^ — 240 nm and A2^^^ 270 nm, re- 
spectively. Taken with a nominal bending modulus of 25 /c^T, 



we estimate the membrane tension to be 4 x 10~^ ksT /nw? . 
From the images, we measure the size of the domains at 
To — 350 — 400 nm, and hence po — 1.5. We estimate the line 
tension, 7, using eqn.fS] based on the fact that the domains are 
dimpled, and find a lower bound of 7 0.49/csT/nm. This is 
in good agreement with theoretical estimates and values deter- 
mined from AFM measurements though somewhat higher 
than the value of 7 ~ 0.22A;sT/nm measured via shape anal- 
ysis of fully phase separated vesicles 0] and 7 ^ 0.40/csT/nm 
from micropipette aspiration experiments [JQ . Finally, viewing 
the repulsive core of the interaction as an effective activation 
barrier to coalescence, a simple Arrhenius argument suggests 
a decrease in coalescence kinetics by two to three orders of 
magnitude. Indeed, such a slowing of coalescence was recently 
observed in a similar model membrane system [44] . 

Discussion 

Comparing biologically relevant domain sizes (^ 50 — 500 nm) 
with the elastic decay length (A2), we expect physiologically 
relevant domains to be small {i.e. small a), as presumed in 
eqn. [H] Estimating the elastic decay length requires knowl- 
edge of the membrane tension and bending stiffness. We 
note that in vitro experiments of osmotically balanced single 
and multicomponent vesicles, and measurements of the plasma 
membrane of unstressed cells suggest membrane tensions of 
10"^ - 10"^ ksT /nm^ gl gD [Ml ES]. The typical bending 
modulus of a phosphocholine bilayer is 10 — 50/cbT, de- 
pending on the exact lipid and cholesterol content [411 1561 [57] . 
Choosing a nominal membrane tension of 10~^/csT/nm^ and 
nominal bending modulus of 25kBT [HJ [J] corresponds to an 
elastic decay length of A2 — 500 nm, suggesting that for lipid 
domains on the order of 50 — 500 nm, small a is an appropriate 
approximation. 

Our experiments on the surface of GUVs have three poten- 
tially confounding effects, all due to the spherical curvature of 
the vesicle. First, the surface metric is not entirely flat with re- 
spect to the image plane. Thus, measurements of distance are 
underestimated the farther they are made from the projected 
vesicle center. This problem is ameliorated by concentrating 
on domains which are at the bottom (or top) of the vesicle 
where the surface is nearly flat and demanding that our track- 
ing software exclude domains that are out of focus; see SI for a 
more detailed explanation. The second potential complication 
is that we use a flat 2D coordinate system for our theoretical 
analysis, however domains reside on a curved surface. Given 
that the domain deformation, and hence energy density, decays 
exponentially with A2 , as long as A2 is small with respect to the 
vesicle radius, the energetics that govern morphology converge 
on an essentially flat surface metric. The final complication 
is that the circular area of focus creates a fictitious confining 
potential for the domains, such that the effective measured po- 
tential of mean force is the sum of the elastic pairwise potential 
and a fictitious potential, Vea — Knt + Vfict • The fictitious po- 
tential is removed by simulating non-interacting particles in a 
circle the same size as the radius of focus (see SI for details). 

The constant tension ensemble used in our theoretical analy- 
sis has a range of validity, determined by the excess area avail- 
able on a thermally fluctuating membrane with conserved vol- 
ume and total surface area Ao (i-e. a vesicle). In the limit 
where the morphological transitions use only a small portion 
(A^) of this excess area, defined by ksT /S-KHih ^ Ao^ the 
tension is constant. Outside this regime the tension rises ex- 
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ponentially with reduction in excess area, tending to stabilize 
dimples from fully budding (see SI for details). 

In addition to the elastic mechanism of interaction, described 
herein, there may be other organizing forces at work in a phase- 
separated membrane: those of elastic ^7], entropic [58l[59] and 
electrostatic origin [60 , however their putative length-scale, on 
the order of ten nanometers or less, is not accessible to the 
spatial and temporal resolution of our experiments, and not 
consistent with our measurement of an interaction length-scale 
of hundreds of nanometers. 

Conclusion 

We have shown that lipid domains are subject to a morpho- 
logical dimpling transition that depends on the bilayer elastic 
properties and domain size. Dimpling allows two domains in 
proximity to repulsively interact due to the deformation in the 
surrounding membrane. Our model makes some key predic- 
tions: at zero spontaneous curvature the domain size distribu- 
tion reaches a critical point where coalescence is arrested by 
repulsive interactions; domains with finite spontaneous curva- 
ture are always subject to interaction and hence should always 
coalesce at a rate slower than the diffusion-limited rate. Ad- 
ditionally, the strength of elastic interactions is augmented by 
increasing line tension or domain area, with an approximately 
linear scaling. The domain size and bilayer elastic parame- 
ters necessary to induce the dimpled morphology are consistent 
with physiological conditions. Further, careful regulation of 
membrane cholesterol in cells may be related to the membrane 
mechanical properties necessary for morphological transitions. 
Combined with lipid recycling, our work offers a mechanism 
working against diffusion-driven coalescence, to maintain fine- 
scale lateral heterogeneity of lipids over time. 

We proposed a simple ID model of an elastic interaction 
that mediates dimpled-domain repulsion, and then used a 
standard ternary membrane system to verify the existence of 
dimpled domains and their subsequent repulsive interaction. 
Finally, it follows that the morphologies and elastic forces 
which organize lipid domains might play an important role 
in the binding and lateral organization of proteins in the 
membrane. 

Materials and Methods 

Giant unilamellar vesicles (GUVs) were prepared from a 
mixture of DOPC (l,2-Dioleoyl-sn-Glycero-3-Phosphocholine), 
DPPC (l,2-Dipalmitoyl-sn-Glycero-3-Phosphocholine) and 
cholesterol (Avanti Polar Lipids) (25:55:20/molar) that 
exhibits liquid- liquid phase coexistence [2 . Fluorescence 
contrast between the two lipid phases is provided by the 
rhodamine head-group labeled lipids: DOPE (1,2-Dioleoyl- 
sn-Glycero-3-Phosphoethanolamine-N- (Lissamine Rhodamine 
B Sulfonyl)) or DPPE (l,2-Dipalmitoyl-sn-Glycero-3- 
Phosphoethanolamine-N- (Lissamine Rhodamine B Sulfonyl)), 
at a molar fraction of ~ 0.005. The leaflet compositions are 
presumed symmetric and hence Vo = 0. 

GUVs were formed via electroformation [2l [61]. Briefly, 
3 — 4 /ig of lipid in chloroform were deposited on an indium-tin 
oxide coated slide and dessicated for ^ 2 hrs to remove excess 
solvent. The film was then hydrated with a 100 mM sucrose 
solution and heated to 50 C to be above the miscibility tran- 
sition temperature. An alternating electric field was applied; 
10 Hz for 120 minutes, 2 Hz for 50 minutes, at ^ 500Volts/m 
over ~ 2 mm. Low membrane tensions were achieved by careful 



osmolar balancing with sucrose ('-^ 100 mM) inside the vesicles, 
and glucose 100 — 108 mM) outside. 

Domains were induced by a temperature quench (see SI) and 
imaged using standard TRITC epi-fluorescence microscopy at 
80x magnification with a cooled (-30 C) CCD camera (Roper 
Scientific, 6.7 x 6.7 /xm^ per pixel, 20 MHz digitization). 
Images were taken from the top or bottom of a GUV where 
the surface metric is approximately flat (see SI). Data sets 
contained ^ 500 — 1500 frames collected at 10-20 Hz with 
a varying number of domains (usually 5 — 10). The frame 
rate was chosen to minimize exposure-time blurring of the 
domains, while allowing sufficiently large diffusive domain 
motion. Software was written to track the position of each 
well-resolved domain and calculate the radial distribution 
function. The raw radial distribution function was corrected 
for the fictitious confining potential of the circular geometry 
(see SI). In the dilute interaction limit, pairwise interactions 
dominate, and the negative natural logarithm of the radial 
distribution function is the interaction potential (potential of 
mean force) plus a constant, as shown in Figure |3]d. 
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